Impact of interactions on human dynamics 
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Queueing theory has been recently proposed as a framework to model the heavy tailed statistics 
of human activity patterns. The main predictions are the existence of a power-law distribution for 
the interevent time of human actions and two decay exponents a — 1 and a = 3/2. Current models 
lack, however, a key aspect of human dynamics, i.e. several tasks require, or are determined by, 
interactions between individuals. Here we introduce a minimal queueing model of human dynamics 
that already takes into account human-human interactions. To achieve large scale simulations we 
obtain a coarse-grained version of the model, allowing us to reach large interevent times and reliable 
scaling exponents estimations. Using this we show that the interevent distribution of interacting 
tasks exhibit the scaling exponents a = 2, 3/2 and a series of numerable values between 3/2 
and 1. This work demonstrates that, within the context of queueing models of human dynamics, 
interactions change the exponent of the power-law distributed interevent times. Beyond the study 
of human dynamics, these results are relevant to systems where the event of interest consists of the 
simultaneous occurrence of two (or more) events. 
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Understanding the timing of human activities is ex- 
tremely important to model human related activities 
such as communication systems [1] and the spreading of 
computer viruses [3]. In the recent years we have expe- 
rienced an increased research activity in this area moti- 
vated by the increased availability of empirical data. We 
now count with measurements of human activities cover- 
ing several individuals and several events per individual 
0i 0i S Hj 01 ■ Thanks to this data we are in a position 
to investigate the laws and patterns of human dynamics 
using a scientific approach. 

Barabasi has taken an important step in this direc- 
tion reconsidering queueing theory 0] as framework to 
model human dynamics [5|. Within this framework, the 
to do list of an individual is modeled as a finite length 
queue with a task selection protocol, such as highest pri- 
ority first. The main predictions are the existence of a 
power law distribution of interevent times Pt ^ t~" and 
two universality classes characterized by exponents a = 1 
(sl-fiol 11 1 and a = 3/2 [l,[lH. These universality classes 
have been corroborated by empirical data for email (5l.[ll| 
and regular mail communications 0, [ill , respectively, 
motivating further theoretical research [1^, |l3|, |lj| . 

The models proposed so far have been limited, how- 
ever, to single individual dynamics. In practice people 
are connected in social networks and several of their ac- 
tivities are not performed independently. This reality 
forces us to model human dynamics in the presence of 
interactions between individuals. Our past experience 
with phase transitions has shown us that interactions 
and their nature are a key factor determining the univer- 
sality classes and their corresponding scaling exponents 
[l5| . Furthermore, beyond the study of human dynam- 
ics, there are several systems where the event of interest 
consists of the simultaneous occurrence of two (or more) 
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FIG. 1: System of two agents with a common interacting task 
I and an aggregate task O representing a set of individual 
tasks. 



events. For example, collective phenomena in disordered 
media, such as the interaction of two (or more) particles 
in cluster formation. 

To investigate the impact of human-human interac- 
tions on the timing of their activities we consider a min- 
imal model consisting of two agents, A and B (Fig. [1]). 
Each agent is modeled by a priority list containing two 
tasks, interacting task (I) and aggregate non-interacting 
task (O). The interacting task models a common activity 
such as meeting each other, requiring the simultaneous 
execution of that task by both agents. On the other 
hand, the non-interacting task represents an aggregate 
meta-activity accounting for all other tasks the agents ex- 
ecute, which do not require an interaction between them. 
To each task we assign random priorities Xij {i = I, O; 
j — A, B) extracted from a probability density function 
(pdf) U,{x) (see Fig.HD. 

The rules governing the dynamics are as follows. Ini- 
tial condition: We start with a random initial condition, 
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FIG. 2: Probability distribution of the interevent time r of 
the interacting task I, as obtained from the direct numerical 
simulations of the model. Each dataset was obtained after 
10^^ model time steps, corresponding with total number of I 
plus O task executions. Note that as La and/or Lb increases 
it becomes computationally harder to have a good estimate of 
Pt because the execution of the I task becomes less frequent. 
The inset shows the distribution for L = 3 as obtained from 
the original model with lO^'^ steps (green diamonds) , and the 
coarse-grained model with A'^ = lO" (red plus), derived to 
obtain more reliable estimation of the exponents. 



assigning a priority to the I and O tasks from their cor- 
responding pdf. Updating step: At each time step, both 
agents select the task with higher priority in their list. 
If (i) both agents select the interacting task then it is 
executed, (ii) otherwise each agent executes the O task, 
representing the execution of any of their non-interacting 
tasks. 

Our aim is to determine the impact of the interaction 
between the agents and the shape of /y [x) on the scal- 
ing exponent a of the interevent time distribution of the 
interacting task I. For simplicity, we focus on the follow- 
ing priority distribution. Consider the case where each 
agent has Lj {j — A, B) tasks, one I task and Lj — 1 
non-interacting tasks, their priorities following a uniform 
distribution in the interval [0,1]. The pdf of the high- 
est priority among Lj — 1 tasks is in this case given by 




FIG. 3: Probability density function of the non-interacting 
aggregate task priority of user A, as obtained from Eq. ([l}. 
With increasing the queue length La, foA{x) concentrates 
more and more in the vicinity of a; = 1~. 



steps. This observation is a consequence of the behavior 
of foj{x) when La and/or Lb are large (Fig. [3]). Focus- 
ing on agent A, as La increases foA{x) gets more con- 
centrated around priority one, while the priority of the 
I task remains uniformly spread between zero and one. 
This fact results in increasingly large interevent times 
between the execution of the I task. 

To speed-off the numerical simulations we derive a 
coarse-grained version of the model, allowing us to an- 
alyze the scaling behavior of the interevent time distri- 
bution over several orders of magnitude (inset of Fig. [2]) . 
We start by noticing that, given {xja, xjb), the joint pdf 
of [xoatXob) factorizes and the probability q{xjA,xjB) 
that both agents execute I right after O is given by 



q{xiA,XiB) 



dxfoA{x) 



dxfosix) . (2) 
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resulting in 



This factorization is possible because the execution of 
the I task requires its priority to be the largest for both 
agents. In turn, with probability 1 — q{xiA, xjb) both 
agents continue to execute O. Thus, the probability dis- 
tribution QrixjA, xjb) that I waits r > 1 steps before 
being executed follows the geometric distribution 
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This example shows that the priorities pdf of task I and 
O are in general different. All the results shown below 
were obtained using the pdf in Eq. 

To investigate the interevent time distribution we per- 
form extensive numerical simulations. Figure [2] shows 
the interevent time distribution as obtained from direct 
simulations of the model introduced above. It becomes 
clear that for large La and/or Lb we do not obtain 
a good statistics, even after waiting for 10"'^"'^ updating 



-{xia,xib) = q{xiA,xiB)[l - q{xiA,xiB)] 
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(3) 



Once the I task is executed it can be executed again re- 
sulting in interevent times of one step (r = 1). The 
overall interevent time distribution of the I task is given 
by 



Pr = 



Pi , T = l 

(1 - Pl){Qr{xiA,XiB)) , T > 1 



(4) 



where 



Pi = 



Si 



(5) 



Si is the expected number of consecutive executions of 
the I task and (• • ■ ) denotes the expectation over different 
realizations of {xja,xib), just at the step of switching 
from task I to O. FinaUy, at the step of switching from 
O to I, the O task priority of both agents must fall below 
that of the I task. Therefore, the pdf of xqj {j — A, B) 
just after the switch from O to I is given by 



foj{x\xio) 



foj{x) 



j;'^ fojix')dx' 



(6) 



where < a; < xjj. This later result together with Eq. 
([3]) allow us to condense all steps with consecutive exe- 
cutions of the O task into a single coarse-grained step. 
More important, this mapping is exact. 

Putting all together the coarse-grained model runs as 
follows. Initial condition: We start with random initial 
priorities extracted from the pdfs fij{x). Updating step: 
At each step, (i) if for both agents the I task priority is 
larger than that for the O task we run the model as de- 
fined above, both agents executing the I task and updat- 
ing their I task priorities using the pdfs fjj {j — A,B). 
(ii) Otherwise, we generate a random interevent time r 
from the probability distribution ([3]) and a new O task 
priority for each agent using the pdf fojix\xij) ©■ This 
second step avoids going over successive executions of the 
O task which, for a large number of non-interacting tasks, 
significantly slow down the simulations. 

The second step of the coarse grained model requires us 
to extract a random number from a geometric distribu- 
tion. This can be achieved very efficiently exploiting the 
fact that the integer part of a real random variable with 
an exponential distribution follows a geometric distribu- 
tion. Using this fact, when t > 1, we extract r exactly 
from the distribution in Eq. which differs from the 
corresponding branch of Eq. Q . Normalization by the 
total number of task I executions, including those with 
T = 1, provides t > 1 distributed according to Eq. 

The I task interevent time distribution obtained from 
simulations of the coarse-grained model is plotted in Fig. 
|4^. When La = Lb = L — 2 it follows a power-law tail 
with exponent a = 2. As L increases a approaches one. 
A guess to this dependence, in good agreement with the 
measured values, is given by a = 1 -t- 1 / max(Lj — 1 ) (inset 
of Fig. 3^). The numerical results indicate that there 
are several numerable universality classes parameterized 
by La and Lb- Notice that the second largest value of 
a (obtained when La = 2 and Lb = 3, or vice- versa) 
is close to 3/2 and, therefore, our results do not show 
universality classes with exponent a between 3/2 and 2 
(unless we assume real valued queue lengths). 
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FIG. 4: a) Probability distribution of the I task interevent 
time for several values of the number of tasks on each queue 
{La, Lb), as obtained from simulations of the coarse-grained 
model. When La = Lb we denote this number by L. The in- 
set shows the exponent a as measured from the power law tails 
(black circles) and the guess function a = 1-1-1/ max(Lj — 1) 
(red curve) in good agreement; to avoid confusion we only 
plot the case when La ~ Lb ~ L, but we checked for the 
general case as well, b) Scaling plot of the I task interevent 
time distribution. Note that, for a given a, the symbols cor- 
responding to different time windows T collapse into a single 
plot. 



The power laws in Fig. exhibit a cutoff at a certain 
value of T. To investigate if this is a natural cutoff or 
just a finite size effect, we investigate the shape of the 
interevent time distribution as a function of the obser- 
vation time window T. The later is defined as the total 
number of steps considering both the I and O task and 
satisfy 



T ■ 



N 
i=l 



(7) 



where N is the number of executions of the I task within 
the time window T and (z = 1, . . . , N) is the sequence 
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of interevent times between executions of the I task. We 
assume that the cutoff is determined by the finite time 
window and that the interevent time distribution follows 
the scaling form 



P(r) = AT-"g 



(8) 



where A is a constant, z > is a scaling exponent and 
g{x) is a scaling function with the asymptotic behaviors 
g{x) « 1 when x ^ 1 and g{x) ^ 1 when x 2> 1. Under 
this assumption P(t) ^ r^" when T oo, with 1 < a < 
2. Given this power law tail and exponent, the number 
of interevent times N necessary to cover the window T 
is of the order of magnitude of T"~^ [l^. In turn, the 
mean intervent time is of the order of 



1 ^ 

{r) = -y 



(9) 



From Eqs. ^ and ^ it follows that z = 1. 

To check our scaling assumption we plot Pt-T" as a 
function of r/T (Fig. |4)d). The symbols corresponding 
to different time windows T clearly overlap into a single 
curve, demonstrating that the scaling assumption in Eq. 
([5]) is correct with z = 1. Thus, in the T — > oo limit the 
I task intervent time distribution exhibits a true power 
law tail Pr ~ r"". 

Within the context of queuing models of human dy- 
namics, only two universality classes were previously 
identified, corresponding to the single queue models of 
Cobham [1,113 (a = 3/2) and Barabasi [1 (a = 1). The 
analysis of the two interacting agents model reveals that 
that the interaction between agents results in a richer set 
of exponents. Our numerical results provide evidence of 
a new universality class with exponent a = 2 and ex- 
ponents between 3/2 and 1. It is worth noticing that 
the exponents 2 and 1 may also result from a Poisson 
proccesses with a time dependent rate Hi 

Because the exponent a depends on the systems de- 
tails, here represented by the agent's queue lengths La 
and Lb, we conclude that the model with two interact- 
ing agents exhibits non-universal behavior. Interestingly, 
the exponent a = 1 is asymptotically reached when the 
number of tasks of one or both agents becomes large. As 
humans get engaged in several tasks this later asymptotic 
behavior may explain the ubiquitous observation of the 



exponent a = 1 11 | 



We use the number of non-interacting tasks as a mean 
to modulate the distribution of the non-interacting ag- 
gregate task priority. Yet, it is the distribution shape 
the primary factor determining the scaling exponent a. 
The effect of increasing La and/or is a concentration 
of the non- interacting aggregate task priority around pri- 
ority one, resulting in values of a that approaches one. 



This means that the limit a = 1 is achieved for low pri- 
ority interacting tasks that remain most of the time in 
the queue without being executed, at expenses of the ex- 
ecution of tasks which in general have a higher priorities. 

Considering the interaction between agents we also 
solve one of the longstanding problems of the original sin- 
gle queue Barabasi model, related to the stationarity of 
the interevent time distribution , . In the Barabasi 
single queue model the task with highest priority is ex- 
ecuted with a probability p, otherwise a task is selected 
at random for execution. When p is close to one the in- 
terevent time distribution exhibits a peak at one step and 
Pi ^ 1 when p ^ 1. When p = I the distribution is non- 
stationary and Pi — > 1 when time t —^ oo. In contrast, in 
the model considered here there is no need to introduce 
the random selection rule and the corresponding model 
parameter p. The interacting task interevent time dis- 
tribution is stationary even when the - highest priority 
first - selection rule is applied. In turn, the exponent a 
is not exactly one, but reaches one asymptotically with 
increasing the number of tasks. Finally, the interevent 
time distribution of the Barabasi model exhibits a natu- 
ral cutoff determined by the parameter p, while for the 
model introduced here it is a true power law up to finite 
size effects. However, it is worth noticing that in a recent 
work fi^l it has been found that the original Barabasi's 
model with variable task execution rate can generate in- 
terevent time distributions with exponent a = 1.25. In 
principle, in our model also different choice of model pa- 
rameters could result in other exponents in the range 
between 1 and 3/2. 

This work represents the first step in understanding 
how interactions among agents affect their activity pat- 
tern. Based on recent works using queueing theory we 
describe the model in the context of human dynamics. It 
can be generalized to consider a larger number of agents 
connected by a specific social network. Also, the model 
can potencially be used more generally to study the time 
statistics of events requiring synchronization between two 
physical systems. 

We thank J.F.F. Mendes for useful comments and sug- 
gestions. 
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